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ABSTRACT 

The satellites of extrasolar planets (exomoons) have been recently proposed as astrobiological tar- 
gets. Since giant planets in the habitable zone are thought to have migrated there, it is possible 
that they may have captured a former terrestrial planet or planetesimal. We therefore attempt to 
model the dynamical evolution of a terrestrial planet captured into orbit around a giant planet in 
the habitable zone of a star. We find that approximately half of loose elliptical orbits result in stable 
circular orbits over timescales of less than a few million years. We also find that those orbits are 
mostly low- inclination, but have no prograde/retrograde preference. In addition, we calculate the 
transit timing and duration variations for the resulting systems, and find that potentially habitable 
Earth-mass exomoons should be detectable. 

Subject headings: planets and satellites: dynamical evolution and stability — celestial mechanics 



1. MOTIVATION 

Exomoons, the satellites of extrasolar planets, have 
been often featured in fiction as habitable locations. 
There is no deficit o f known giant planets; Exoplanet.org 
(| Wright et al.1 120101 ) lists approximately 40 giant exo- 
planets (8% of total) within 20% of the equilibrium tem- 
perature of Earth, as are 30 (3%) of the Kepler planet 
cand idates released in February 2011 (|Borucki et al.l 
2011). Though these observations are preliminary, they 
do show that habitable-zone giant planets not only ex- 
ist, but are common. Once a giant planet is known to 
be in a habitable zone, variatio ns in its orbit, such as 
Tran sit Timing Variation (TTV; ISartoretti fc Schneider 
1999) and Transit Duration Variatio n (TDV; IKippme 
2009), photo metry dSzabo et al.1 120061). or gravitational 
microlensing (|Liebig fc Wambsganss 2010]), allow the in- 
direct detection of satellites. Thus, if potentially hab- 
itable exomoons exist around transiting giant planets, 
they may be detected at the same (or even greater) rate 
as solitary habitable terrestrial planets. As yet, no exo- 
moons have been detected, but the wealth of transit data 
from the Kepler mission should begin to fill this gap. 

Despite the existence of giant planets in stellar habit- 
able zones, it is far from certain how they arrived there. 
Current giant-planet formation models assume that they 
are created at distances beyon d the stability point of ice 
(e.g. ILissauerlll987t Boss 1997|), which implies conditions 
not suitable to surface habitability. D isk migratio n can 
bring giant planets close to the star (|Wardl Il9~97j) . but 
generally has a stopping point far too close to the star 
to be habitable (thus producing "Hot Jupiters"). The 
host planets of potentially habitable exomoons there- 
fore likely arrived at their final orbit throug h late-stage 
migra tion, driven either by planetesimals (|Kirsh et aL 
2009) or other giant planets (jW eidcnschil ling &: Marzari 
1996). 
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In the process of migrating, the satellite systems of 
these giant planet may have close encounters with ter- 
restrial planets or planetesimals, causing them to be dis- 
rupted or replaced. If either the Jovian or Saturnian 
systems were transported to 1 AU around a solar mass- 
star, both Callisto and Titan would be at 18% of their 
planet's Hill radii, thus implying that all the major satel- 
lites of the two planets would be on stable orbits. How- 
ever, a close encounter could either excite their orbits to 
high eccentricity (thus requiring tidal recircularization), 
or could result in the capture of a much larger terrestrial 
satellite. Neptune appears have to experienced this pro- 
cess during its migration through the proto-Kuiper Belt, 
loosing any original major satellites, while gaining Tri- 
ton in an inclined, retrograde orbit. This was possibly 
due to a momentum-exchang e reaction that ejected th e 
binary companion of Triton (jAgnor fc Hamilton! F2006) . 
though other scenarios are possible (at reduced probabil- 
ity). Any capture process, though, will tend to produce 
very loosely-bound initial orbits, with only a small delta- 
v to escape velocity at periapse. Therefore, some method 
must be used to determine the long-term evolution and 
stability (or lack thereof) for these orbits. Here we use 
a full KCTF (Kozai Cycle and Tidal Friction) model to 
find the survival probability for a range of physical con- 
ditions and the d etectability of the result ing system. 

As shown by IDonniso n| (pOlOl ) and ISato k Asadal 
(2010), there are a range of stable orbits for Earth-mass 
planets around giant planets. However, both of those 
models test only the stability of the orbit, rather than 
any evolution due to tidal effects. On inclined exomoon 
orbits, the effects of stellar torques on the orbit can, 
through initiating Kozai cycles, dramatically accelerate 
the rate of tidal decay, orbit circularization, and spin- 
orbit synchronization. As we will show, this process can 
allow even very loose, inclined capture orbits to stabilize 
to tight, circular orbits. Thus, even marginal dynamical 
captures of former terrestrial planets can produce stable 
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Fig. 1. — Post-capture spin-orbit evolution of two exomoons; the 
solid line shows a Kozai-enhanced decay, while the dashed line 
shows a non-Kozai decay. 



exomoons. 



2. KCTF MODEL 

In order to understand how exomoon orbits may evolve 
after being captured, we created a numerical Kozai Cy- 
cle and Tidal Friction (KCTF) model. Kozai cycles in 
this context are the secular oscillations in eccentricity 
and incl ination of th e exomoon's orbit caused by stellar 
torques (|Kozailll962h . In isolation, these oscillations pre- 
serve the semimajor axis of the orbit and the quantity 
Hk = cos(I) x Vl — e 2 , where I is the inclination of the 
exomoon's orbit relative to the planet's stellarcentric or- 
bit, and e is the eccentricity of the exomoon's orbit. In 
the scenarios simulated, these cycles can initially be as 
fast as just a few years (see Figured]), causing very rapid 
evolution of the exomoon's orbit. Since all the initial or- 
bits for this study were highly elliptical, only very low 
initial inclinations (within about 15° of coplanar) pro- 
duced values of Hk sufficiently low that Kozai cycles did 
not occur. 

Because the Kozai oscillations from the star attempt 
to preserve Hk, a drop in the inclination of the orbit 
can cause the eccentricity to become very high. As the 
eccentricity of the orbit increases, the periapse of the ex- 
omoon's orbit becomes much closer to the planet. Tidal 
forces become much stronger at these close distances (go- 
ing as a~ 8 in this model), and so a close periapse due 
to eccentricity causes the orbit to shrink and thus decay 
further. This sets up a positive feedback loop, which pro- 
gressively shrinks and circularizes the orbit. The stair- 
step semimajor axis decay in Figure Q] happens because 
the obit is initially only decaying at high eccentricities; 
once the apoapse is close enough to also experience tidal 
forces, the oscillations stop. The result is to cause semi- 
major axis decay and circularization much faster than 
if the star were not exciting eccentricity. In addition, 
since the high eccentricity part of the Kozai cycles are at 
low inclinations, orbits are preferentially frozen near the 
plane of the stellarcentric orbit. 

To simulate th is process, we used a KCTF mode l 
based on that of lEggleton fe Kiseleva-Eggletonl (|2001[) . 
which allows the direct integration of the exomoon or- 
bit's specific angular momentum vector h and Laplace- 
Runge-Lenz eccentricity vector e, as well as the spin vec- 
tors of the planet and satellite. The tidal properties 
of the giant planets were based on those presented in 



iFabrvckv fc Tremaind (|2007f) for a Jupiter-mass planet. 
The tida l properties of th e exomoons used the formu- 
lation of Ragozzinc (2009), along with his addition of 
solid-body quadrapole gravity. To integrate the system 
of equations, we used the B urlisch-Stoer met hod with 
vector-rational interpolation (Sweatman 199 a) and er- 
ror co ntrol based on the algorithm given in I Press et al.l 
(2007]), with a per-timestep tolerance of 10~ 10 . 

This model does not include any dynamical effects from 
objects external to the exomoon-exoplanet system other 
than the star itself. In addition, we assumed that the 
planet did not migrate over the course of the simulation. 
If it were to migrate, the effect would be to shrink the Hill 
radius of the planet and shorten the period of the Kozai 
cycles. The smaller Hill radii would allow more captured 
satellites to escape, while on the other hand, the faster 
Kozai cycles would decrease the decay timescale. Thus, 
if this model works fast enough for a static case, it should 
also be applicable to a slowly migrating planet. 

3. MONTE CARLO SIMULATIONS 

To find in what conditions a captured exomoon may 
survive, we generated 18 sets of synthetic Star-Planet- 
Moon systems and performed KCTF simulations on 
them. Each set contained 1000 synthetic systems, with 
common masses for all objects, and randomized initial 
exomoon orbits and spin states. To simulate a loose cap- 
ture, the initial exomoon orbits all had apoapses beyond 
80% of the planet's Hill radius, and eccentricities greater 
than 0.85. This is most consist ent with a low mass ratio 
momentum-exchange capture (jFunato et al.l I20041 ). but 
generally similar to any low delta-v, non-disruptive cap- 
ture (or eccentricity excitation due to a close encounter). 
Again consistent with a capture, both the satellite's ini- 
tial orbital plane and spin vector were initially pointed 
at random directions on the sky. This means an approxi- 
mate equipartition between prograde and retrograde ini- 
tial orbits and between prograde and retrograde initial 
spin states for the satellite. The planets had a random 
obliquity less than 5° from their stellarcentric orbit. Each 
planet-moon system was at a stellarcentric distance such 
that the equilibrium temperature was equal to Earth. 
The stars and stellarcentric orbits used were the Sun 
(G2) at 1.0 AU, a main-sequence F0 (1.7 M Sun ) at 2.1 
AU, and a main-sequence M0 (0.47 Msun) at 0.28 AU. 
The planets we assumed to have a mass equal to either 
Jupiter or Neptune, u s ing th e tidal parameters given in 
IFabrvckv fc Trema ine (2007), though nearly all dissipa- 
tion was in the exomoon. The simulated exomoons had 
the mass of either Earth, Mars, or Titan (with Mars un- 
compressed density), usin g a tidal Q of 100, mo dulus of 
rigidity of 3 x 10 10 7V/m 2 (jCladman et al.lfl99"rl . and J 2 
of 0.001. The initial rotational periods of the planet and 
satellite were varied randomly between 0.1 and 48 hours. 
The simulations were run until they reached either 10 9 
years or an eccentricity below 10~ 5 . However, the simu- 
lations were terminated early if the periapse went below 
the Roche limit (thus potentially causing breakup of the 
exomoon) or the apoapse exceeded the Hill radius (allow- 
ing the exomoon to become unbound) . It is possible that 
the exomoon could have survived and remained bound to 
the planet in these scenarios, but that is beyond the fi- 
delity of these simulations. 

To estimate the detectability of the resulting systems, 
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FlG. 2. — Post-evolution orbital period and inclination distribu- 
tion; most orbits are coplanar with the planet's orbit, with no 
pro/retrograde preference. 
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FlG. 3. — Post-capture circularization timescalcs; note that nearly 
all are below a few million years, and timescales decrease with 
stellar mass. 



we calculated root-mean-squared Transit Timing Varia- 
tions (TTVs) and Duration Variations (TDVs) for each 
simulation. The TTV and TDV are due to the orbit of 
the exomoon causing the exoplanet to begin the transit 
either sooner or later than the barycenter of the planet- 
moon system. These effects are maximized for low mass 
ratios (e.g. Earth mass exomoon around a Neptune- 
mass planet) and longer period exomoon orbits. Since 
the majority of resulting systems were low-inclination 
with respect to the stellarce ntric orb it , we u sed the zero- 
inclination equations from iKippinj (|2009( ). Assuming 
zero-inclination maximizes the TTV and TDV, allowing 
an estimate of whether the system would be even de- 
tectable in the best case scenario. 

4. RESULTS AND DISCUSSION 



The general result of the simulations is that loosely- 
captured exomoons around giant planets in habitable 
zones can survive to stabilize into long-lived orbits. Table 
Q] shows the fraction that stabilize for each mass/distance 
scenario. The exomoons stabilized much easier in the 
FO case than the MO case, with the solar-mass case in- 
between. The reason for this is readily apparent from the 
observation that while Hill radius scales linearly with dis- 
tance from the star, the equilibrium temperature scales 
as the inverse square of distance. Thus, the super-solar 
mass star allows a much wider Hill radius at a habitable 
distance than for a planet of equal mass around a sub- 
solar mass star. This larger Hill radius provides the post- 
capture orbit much more room to move around, allowing 
longer period exomoons to stabilize. The truncation of 
periods greater than two days for the M0 case in Figure 
[3] is a clear result of this. 

However, it also apparent from Figure [3] that the fast 
stellarcentric orbits of the MO case allow for much more 
rapid exomoon circularization that for the larger stars. 
This clearly shows the effect of Kozai-pumped eccentric- 
ities accelerating tidal decay, as otherwise the timescales 
would be independent of the distance from the star. The 
median circularization timescales for the M0 cases are of 
order 10 4 years, 10 5 years for the solar-mass star, and 
10 6 years for the FO star. All of these timescales are very 
short relative to the lifetime of the star, but 10 6 years 
may be long enough that the exomoon's orbit could be 
externally perturbed by a planet or planetesimal disk. 
On the other hand, as Figure Q] shows, the semimajor 
axis decay typically happens at least an order of mag- 
nitude faster than full circularization. Therefore, it is a 
reasonable assumption that nearly all these orbits would 
stabilize before any external perturbation would disrupt 
them. 

To gauge the amount of tidal heating on the satel- 
lite, we estimated the amount of energy lost from the 
system from the difference in the initial and final orbits 
and spins. For the orbits, we first found the energy dif- 
ference between the initial and final states. Then, we 
muliplied this an estimate of the fraction of energy that 
went into the satellite using the masses and tidal Q of 
the objects, Q p ianetM sat /(Q p i anet M sat + Q sa tM p i anet ). 
This term was usually near unity, as the planet was as- 
sumed to have a very large Q. We then estimated the 
change in rotational energy of the satellite, subtracted 
this from the change in orbital energy, and divided the 
result by the circularization time to produce a heating 
rate. The rate was dominated by the orbital term, as 
the initial rotation periods were not too different from 
the final orbital periods. Generally, this rate was higher 
per unit mass for larger satellites (which circularized 
faster) and closer orbits. For the Earth-mass/ Jupiter- 
mass case, the median rates were 0.002 mW/kg around 
a FO star, 0.02 mW/kg around a solar-mass star, and 0.5 
mW/kg around a M0 star. These rates are higher than 
for either long-lived (w 10~ 8 mW/kg: iDesch et alJl2"009[ ) 
or short-lived chond ritic radioisotopes (~ 10~ 4 mW/kg; 
ISchubert et H1I2007I ). A more through analysis is be- 
yond the scope of this work, but these rates could create 
significant short-term heating on the exomoons. 

Table [T] also shows the fates of the exomoon orbits 
which did not survive to circularize. The systems around 
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TABLE 1 

Relative fraction of end states for fully evolved exomoon systems. 



Star 


Planet 


Moon 


Survived 


Retrograde 


Separated 


Impacted 






Earth 


43% 


52% 


21% 


35% 




Jupiter 


Mars 


44% 


45% 


18% 


37% 


Sun 




Titan 


42% 


47% 


21% 


36% 




Earth 


52% 


44% 


17% 


30% 




Neptune 


Mars 


44% 


45% 


18% 


36% 






Titan 


45% 


47% 


19% 


35% 






Earth 


65% 


47% 


3% 


31% 




Jupiter 


Mars 


59% 


46% 


4% 


35% 


FO 




Titan 


61% 


48% 


3% 


34% 




Earth 


77% 


44% 


4% 


18% 




Neptune 


Mars 


67% 


44% 


4% 


28% 






Titan 


61% 


50% 


4% 


33% 






Earth 


23% 


50% 


2% 


74% 




Jupiter 


Mars 


23% 


51% 


3% 


73% 


MO 




Titan 


23% 


44% 


2% 


74% 




Earth 


24% 


50% 


1% 


73% 




Neptune 


Mars 


25% 


52% 


1% 


73% 






Titan 


22% 


50% 


3% 


74% 



TABLE 2 

Median Exomoon orbital period, Transit Timing Variation (TTV), and Transit Duration Variation (TDV) for full-evolved 

exomoon systems. 



Star 


Planet 


Moon 


Period (d) 


TTV (min) 


TDV 


edge-on (min) 






Earth 


2.17 


5.44 


0.114% 


0.93 




Jupiter 


Mars 


2.34 


0.59 


0.012% 


0.10 


Sun 




Titan 


2.52 


0.12 


0.002% 


0.02 




Earth 


1.65 


36.69 


0.835% 


6.41 




Neptune 


Mars 


1.89 


4.11 


0.089% 


0.69 






Titan 


2.16 


0.86 


0.018% 


0.14 






Earth 


3.68 


11.33 


0.111% 


1.45 




Jupiter 


Mars 


4.26 


1.22 


0.011% 


0.15 


P0 




Titan 


4.14 


0.26 


0.002% 


0.03 




Earth 


2.38 


75.70 


0.860% 


10.60 




Neptune 


Mars 


3.79 


8.51 


0.083% 


1.03 






Titan 


3.50 


1.78 


0.018% 


0.22 






Earth 


0.94 


1.63 


0.122% 


0.54 




Jupiter 


Mars 


1.00 


0.18 


0.013% 


0.06 


M0 




Titan 


1.00 


0.04 


0.003% 


0.01 




Earth 


0.85 


10.97 


0.836% 


3.36 




Neptune 


Mars 


0.79 


1.22 


0.098% 


0.39 






Titan 


0.83 


0.26 


0.020% 


0.08 



the solar mass star were by far the most prone to separat- 
ing beyond the Hill radius, with roughly equal impacts 
as separations. The M0 cases, on the other hand, actu- 
ally impacted the planet in three-quarters of the cases, 
due to the very small Hill radius not allowing the pe- 
riapse to shrink very far before impacting. A majority 
of F0 cases survived in all scenarios, with the large Hill 
radius allowing most orbits to circularize. Therefore, for 
a given capture rate, potentially habitable captured ex- 
omoons should be the most common around stars larger 
than a solar mass. 

Figure [2] shows that the majority of the resulting or- 
bits, especially around a Jupiter-mass planet, are nearly 
coplanar with the stellarcentric orbit. This is partially 
due to the assumed low obliquity of the planet align- 
ing its gravitational quadrapole with the direction of 
the star. However, as shown by the high inclination of 
Triton (5.9 day period), an orbit does not need to be 
very large for this effect to be minimized. And indeed, 
the orbits around Neptune-mass planets show a much 



larger diversity of inclinations for periods less than two 
days, especially for Earth-mass satellites. Beyond this 
point, solar torques begin to dominate and force nearly 
all orbits coplanar. Since longer period orbits produce 
larger TTVs, it is therefore likely that the first exomoons 
detected will be coplanar with the stellarcentric orbit, 
and it will take much more observations to detect short- 
period, high-inclination exomoons. 

Also congruent with Triton, Figure [5] shows that there 
is no real preference for either prograde or retrograde 
orbits. The secular perturbations of Kozai cycles do 
not allow the directionality of the orbit to change in 
all but the most inclined initial cases. Thus, each or- 
bit preserves its original direction, evolving towards the 
plane of the stellarcentric orbit. All of the final orbits 
are circular and have semimajor axes less than 28% of 
the pla net-moon system 's Hill radius, which as summa- 
rized bv lDonniso means that they shoul d be sta- 
ble oy er the lifetime of the star. Unfortunately, [Kipping 
(2009) shows that TTV and TDV are both degenerate 
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for prograde/retrograde inclinations, and so measuring 
the relative fractio ns may be very difficult. However, 
iSimon et all (|20f0l ) show that it may be possible to break 
this degeneracy with very precise measurements of the 
Rossiter-McLaughlin effect. Since secular perturbations 
cannot reverse the direction of a satellite's orbit, detec- 
tion of a retrograde satellite around an exoplanet would 
be confirmation that this capture process occurred. 

The median TTV and TDV for each simulation set 
are given in Table [T] The TTV do vary over a large 
range, but the upper end of that range is just within 
the limits of detection for modern transit systems. The 
Kepler mission has a regular cadence of approximately 
every 30 mi nutes, and can run at a cadence of up to once 
per minute (jKoch et alj|2010j ). This appears sufficient to 
detect an Earth-mass exomoon in most cases, and just 
enough to detect a Mars-mass exomoon in the best cases. 
Larger stars offer wider orbits (and thus stronger TTVs), 
but much more infrequent transits. Solar mass to early 
K stars therefore seem to offer a good balance between 
TTV strength and transit frequency. 

The T r ansit Duration Variation (TDV) introduced by 
iKippingj (|2009f ) allows the unambiguous detection of an 
exomoon from transits alone. For the systems described 
in that work, the TDV was of similar duration to the 
TTV. However, we find that the lower planetary masses 
and much wider stellarcentric orbits of our study result 
in relatively weak TDVs for nearly all cases. Thus, it 
is unlikely that a potentially habitable exomoon will be 
detected by TDV alone; TDV could, though, allow for 
follow-up transit observations to confirm the existence of 
a TTV-detected exomoon. 

These results compare well with previous estimates, 
though both of those simulations assumed much wider 
orbits than were usually produced by the KCTF model. 
The consequence of this is both reduce the TDV sig- 
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nal relative to that a ssumed bvlKippind (2009J) and the 
photometric signal of lSzabo et al.l (|2006f) . On the other 
hand, Table [T] shows that even with these closer orbits, 
some exomoons are still within the range of detectability. 
The combination of TTV and TDV may offer a stronger 
detection signal than photometry for these orbits, though 
both could detect some of the orbits produced. 

5. CONCLUSIONS 

Several general conclusions can therefore be drawn 
from our simulations: 

f. Loosely-captured exomoons around giant planets 
in habitable zones can survive to stabilize into long- 
lived orbits. 

2. The timescale for them to stabilize is enhanced by 
Kozai torques, and is generally less than a few mil- 
lion years. 

3. Most of the surviving orbits are close to the 
plane of the exoplanet 's orbit, but show no pro- 
grade/retrograde preference. 

4. Some of the resulting orbits should be currently de- 
tectable (especially for Earth- mass satellites), with 
the transit timing variation much stronger than the 
duration variation. 

5. The most promising targets for detecting poten- 
tially habitable exomoons by TTV appear to be 
Neptune-mass exoplanets orbiting stars of solar 
mass or slightly below. 
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